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Abstract 

The concept of the temperature of a surface is introduced from the viewpoint of the 
physical chemistry of surfaces. The surface, near surface and microlayer regions of the 
interface are defined. Most methods measure the temperature of the microlayer or at best 
the near surface region and may err in representing the surface temperature. Methods based 
on capillary ripples actually measure the surface temperature since surface tension (or 
surface tension tensor when a monolayer has been spread or adsorbed at the interface) is the 
main restoring force that controls their propagation. Light scattering methods are described 
for determining the elevation of very small amplitude capillary waves through the 
computation of various correlation functions from which the surface tension can be 
estimated. Procedures for estimating the surface temperature are described. 


The Surface Temperature 

In order to speak of the "temperature of a surface" a length scale for the surface must 
be specified: the interface spans a length over which the normal derivative of the density is 
non-zero, Figure 1. It is well understood theoretically and experimentally that away from 
critical temperatures, the width of the interface is at most a fraction of a nanometer and 
corresponds to a few atomic or molecular dimensions. (See ref. (1,2).) 

It is useful to define two additional length scales. The near— surface extends to about 
10 nm. The near surface properties are some superposition of surface and volume properties 
which may, with careful analysis, be separated to some extent. The microlayer can be 
thought of as extending to roughly 1 micrometer from the surface. The properties 
determined in the microlayer cannot be separated into distinct surface and volume effects. 

A property such as the surface temperature refers only to the first few molecular 
(atom) diameters of the interface; rigorously, it is not either the near-surface temperature or 
the microlayer temperature. Most methods in fact report the near surface or microlayer 
temperature. When the system (e.g. two-phase with one component) is at equilibrium, the 
temperature is uniform so that T' = T gurf = T"; even the microlayer temperature 

represents the temperature of the entire system. When the system is not in equilibrium, the 
near surface or microlayer temperatures may differ significantly from the surface 
temperature. The hypothesis that the near surface and surface temperatures are the same 
must be established experimentally or supported by careful analysis. 
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Ellipsometry, reflectrometry (especially in the x-ray wave lengths) as well as theory 
and computer simulation provide the basis for asserting that the liquid/vapor interface is at 
most a few tenths of nanometers in thickness. One representation of the density profile 
through the interface is that of the mean field, Figure ( 1 ): 
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where n is the volume density, n^ + \ n^ ^ the densities on either side of the interface, and z is 
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value of £ depends on T. In addition the interface fluctuates. Define the Gibbs surface 
associated with the interface by the integral 
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where -a and a are taken so that (n(± a) -n*) ~ 0 . 


The fluctuating interface can be represented by the differential manifold 
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One must think of ( as a random variable, for the plane interface, or as a small fluctuation 

2 

away from the average trace of a curved interface. The representation for a plane interface 
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for unit area leads to 
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(for unit area of surface) where k is the Boltzmann factor, T the temperature, 7 the surface 

tension, q the wave number, An = n^ - n^ ^ the density difference and g the gravity 
constant. This formula has been studied experimentally and is considered accurate although 
a molecular level theory is incomplete. 


* 

Clearly, the measurement of will provide an estimate of the true surface 

temperature. However, the surface tension, 7 , depends on temperature. So that the 
temperature dependence of the surface tension should be known. 

In addition, the spectrum of the fluctuations can be calculated (Ref ( 2 ) and references 
therein) 
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( 6 ) 


G(w,q) = F(< Cg(0) C-(t)>) 

* -> 
where <C^(0) C^(t)> is the autocorrelation function of the q component of the elevation 

function and F( • ) is the Fourier transform. 


Even though J <( 2 > << -1 nm, continuum mechanics provides a formula for 

which is known to be accurate for small |q| (the "long" wave length limit). There is a 


substantial body of experimental work that supports this conclusion so long as |q| is 
sufficiently small that macroscopic formulas are valid. A fairly delicate argument leads to 
the following result (Ref (2) and references therein) 


G^{u) = - % kT ( M bii [q,^,n,//,7,. . .] (7) 

4 7q + gAn ' 

where , . . ., represents surface visco-elastic coefficients appropriate for systems for which a 
monomolecular film can form; a one component, two phase system admits only 7. 

Note that the determination of G^(w) provides the surface temperature as well as the 

surface tension. However, there are conditions that must hold when the system is dynamic, 
and VT / 0. The distribution function for {<%} must not change form from that of 

equilibrium when the system is dynamic; local equilibrium is assumed. Furthermore a 
sufficient area must be sampled so that area averaging is satisfactory; for an t by l surface 

patch, (» i where q ~ 1000 cm~^. 


* 

In fact, the functions <(qCq> and G (w) can be determined by light scattering which 
will be described briefly. 


Light Scattering Spectroscopy 

Early in the century, the Maxwell equations for the scattering of light from thermally 

excited capillary waves 2,3 was worked out with the result that in the far field for C-» « A , 

b q o’ 

with A q the wave length of the incident light, the scalar electric field follows E-> = (% x f 

where f depends on the geometry of the scattering experiment (angles of incidence and 
scattering) and the refractive indices of the phases that join at the surface. Obviously, at 

least one of the phases must be transparent. Indeed <E^ E^> = r<(-> (-*>. However, the 

electric field is detected by a square-law device (photodiode or PMT detectors) so that the 

autocorrelation function actually generated is that of the current, i = E E . Thus a 

heterodyne system is required so that <i(0) i(t)> can be related directly to <<^(0) C _+ (t)> 

o a c q q 

without loss of phase information 
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A general schematic of the methodology ’ ’ is shown in figure 2. The essential 
feature is that a spatially filtered image of a grating of known wave number is projected onto 
the interface. The quality of the optical processing in projecting the grating onto the surface 
is crucial. The liquid surface is usually the highest quality optical component in the system 
and can be thought of as a very high quality but very inefficient phase grating following the 

rule 4 ' 5 

R = R 0 (l + (8) 


where R is the reflection coefficient of the surface, R is the reflection coefficient of the 

ideally smooth surface, r/(x, y) << 1 is the amplitude variation over the illuminated spot and 
<£{x,y) is the phase variation. The phase function is <^j(x,y) = 2k Q <^+ wherein a real 


function representation of is appropriate 


4,5 


The product RS^ ^ provides the modulation of the reflected beam, by the 

capillary ripples; a particular component (j^(t) modulates (= RS^ depending on the 

structure of ^ as determined by the grating, F^, and L 9 . However, since the incident 

laser beam is of finite extent, the reflected beam carries a distorted representation of 

(^(t). Fortunately, the "instrument function" can be computed accurately^. The one most 

characteristic dimensionless group for the instrument function is k a where k is the wave 

g g 

number of the grating as projected on the surface and a the cross section of the Gaussian 
incident beam (4 a is the "beam diameter"). While not absolutely essential, the incident 
beam should be cleaned up so as to be accurately Gaussian since this makes the computation 
of the instrument function relatively simple. 


It is practical to use optical components to manipulate the optical field so as to 
eliminate the zero-order beam and pass only the reference beam that is modulated by <p. 
The optical quality required for this step is modest, the entire reference beam must be 
collected excluding all other beams (higher order or lower order beams). The beam need not 
be of accurate Gaussian cross section and some distortion of the phase by the collection 
optics can be tolerated. 


The computation of the correlation function <i(0) i(t)> = R- (t) from the time series 

* 4 5 

generated by the detector yields <E^(0) E^(t)> and therefore ’ 



dq y <C*(0) C q (t)> F(q x -k“,q y ) 
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where k 


s 

g 


is the wave number of the grating as projected on to the surface (if the 


magnification is 1 then k 


g 


c c 

k ). The instrument function F depends on k a so that 
g g 
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k g a— 4 °°’ F “ 4 ^ q x ” k g) S (Qy) and R j — 1 < < q (0) < q (t)>, see references (4,5). 

A second dimensionless group useful for designing an experiment^’^ is Y = (u 0 t 0 )^ 

where = 2-- j^, + where fi'fi ' ' are the shear viscosities in the volume of the 

phases surrounding the interface and n', n' ' the densities. The frequency, u 0 is defined by 
z 2 3 

1 = (n' + n' / )a; Q /7q which obtains for the ideal, clean surface, between fluids of 
vanishing viscosity. 

When Y >> 1 the capillary ripples are propagating and 

<C(0) C(t)> = e -r q l cos w t (11) 

q q 7Q + gAn q 

is a good representation of the autocorrelation function. The correlogram determined 
experimentally allows the estimation of and (and their uncertainties and correlation) 

by fitting equation (11) through a maximum likelihood algorithm. In general is 

determined to ~ .2% and to ~ 1%. It then is easy to compute any two parameters n,/j,7, . . 

. by using the dispersion equation, Det M = 0, from equation (7). 

The construction of equation (7) and in particular the response function matrix M, has 

been done several times and is too lengthy to produce here; see references (2,3) for a guide to 
the literature. The analysis requires a careful construction of the mass and momentum 
balance at the surface including coupling to the substrate through the jump of the three 
dimensional pressure tensor across the interface projected along the normal to the interface 

o o 

as well as in the plane of the interface ’ . The boundary conditions include: 



where v | Q is the volume velocity field evaluated at the surface and projected either along the 

normal, A, or into the plane by the tensor U The conditions for linearization obtain since 

C << ^ i n the spectrometer specifications. M is constructed from the set of linear equations 

that result from a Fourier trahsform in space and a Laplace transform in time. The M7j 
component representing ( is used in equation (7) while Det M = 0 is the dispersion equation 
mentioned above. 
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Remarks on Implementation of the Fourier Transform Spectrometer 

A research instrument, as figure (2) suggests, has required substantial room on an 
optical bench. Optics of focal lengths of ~ 1 meter are often used. We have built a much 
more compact system for various applications. Our latest instrument uses a ~ 40 mw diode 
laser, A ~ 830 nm, which is only a few centimeters long, wide and high. The power supply is 

built into this package. With care, a lens system of f - 200 mm works well to provide the 
necessary Fourier transforms and further, the optical path can be folded. The detector is a 
photodiode and is small. This system delivers correlograms of research quality without 
special manipulation of the output of the laser; the beam is elliptical but of Gaussian 
cross section in each direction. We are now working to reduce the size to that of a package 
measuring centimeters instead of meters. 

A second problem with the light scattering methodology has been sensitivity to 
vibration. When the surface tilts slightly at low frequency (< 10 Hz ) due t0 building 

vibration the correlograms become so badly distorted that T , cannot be determined. 

This has been overcome by a beam directing device in a feed-back circuit that locks the 
reference beam onto the detector accurately. The correlogram is then clean of the flicker 
effect; the slight distortion of the correlogram due to fluctuations of the angle of incidence ol 
the laser beam about the normal is ignorable. 

We have designed and are implementing a "single board" correlator that is driven by a 
PC (80286 processor + 80287 coprocessor) which provides a compact, clean system for data 
analysis. 

We expect to have the first integrated prototype ready for a ship-board experiment in 
October, 1989 involving the measurement of the calm ocean surface tension in situ. 


Discussion 


The light scattering procedures for measuring the "true" surface temperature rather 
than the near-surface or microlayer temperature are two 


• Measure the averaged, surface light scattering irradiance which is proportional 
to <C ( > and thereby determine T through equation (5). 

• Measure the correlation function of the scattered light, compute estimates of T 
and uiq and determine T through their variation with T. 


The determination of the averaged, surface light scattering irradiance is difficult 
because the signal is small and subject to distortion by flare. Moreover, even with the 
careful application of an interference filter with a narrow band-pass at the wavelength of the 
laser, the radiance of the heated, liquid metal pool may swamp the detector. However, since 
the laser can be modulated easily, lock-in detection is possible and should be effective in 
separating out the signal due to scattered light. The strong q dependence of the scattered 
irradiance also provides a way of separating out the part of the signal to be analyzed. These 
methods require some information about the equilibrium surface tension dependence with 
temperature but a limited data set should be sufficient because of the relatively small 
changes involved. 
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The measurement of the correlation function is quick and accurate. The radiance 
coming from the heated, liquid metal pool is uncorrelated but does contribute to the "dark" 
noise. An interference filter should make this effect tractable. Note that the reference beam 
is much more intense than the scattered light and so the signal-to-noise ratio will be 
substantially better than in the first procedure. 


However, the relationship between T , and T is complex in that the surface tension 

dependence on temperature is involved as well as that of the density and viscosity. Even so, 
the methodology has the potential of being relatively easy to calibrate (w , with T for 


each material) under equilibrium conditions which provide working curves for determining T 
and 7 for non-equilibrium systems. 


While the theory of light scattering spectroscopy can be generalized to include 
temperature gradients, at this time we depend on the local equilibrium assumption. If light 
scattering can be observed from a small enough spot then local equilibrium can be asserted 
and an estimate of the surface tension or temperature, averaged over the illuminated spot, 
can be computed. However, as the spot is made smaller, q becomes less well determined by 
the optical system and thereby the estimate of the surface tension and temperature is less 
certain. 


With the technology we have in place now, the minimum spot size that is practical is 

about 2 mm. It is likely that we can work with .2 mm spot sizes by using gratings with k > 

§ 


1000 cm 1 but how well this works in practice must be determined. Note: The spot can be 
moved over the surface so that a profile can be developed. 


It is entirely possible to do surface spectroscopy from curved surfaces. Indeed little 
needs to be changed for cylinder systems of 1 cm or so in length and a few mm in diameter or 
for spheres of 1 cm or so in diameter. It is an interesting question whether by careful 
application of Fourier transform optics it is possible to design an optical system that will 
separate out the higher frequency normal modes of a spherical drop for spectral analysis. 
Based on unpublished calculations done for spherical drops irradiated by an ultra-sound field 
we feel it is possible. 


The question is whether an optical field can be produced in practice so that normal 
mode oscillation driven by thermal fluctuations in a small drop can be detected as a time 
series and analyzed. Indeed, it is not necessary to consider thermal fluctuations alone, 
external fields (e.g. sound or ultra-sound) can be used to drive the surface oscillations of the 
drop. The spectroscopy of such driven waves is also possible using the laser methods 
developed by the CWRU group. 

It has been our experience that the light scattering methodology is the method of 
choice in determining the surface tension and interfacial tension of a system. It is fast, 
accurate and does not require manipulation of the sample. It is my opinion that this is the 
method of choice for determining the equilibrium surface tension of melts at high 
temperature. Whether the technology can be adapted for small samples with curvature can 
be determined by Fourier optical computations and a set of careful measurements. 
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Figure 1 The density variation through the interfacial region between a liquid and its 
vapor, a) The variation of density as determined by theory and simulation, b) 
A model that can be defined to be consistent with the macroscopic continuum 
thermodynamics of interfacial behavior. In either case, the interfacial thickness, 
£, is < 1.0 nm for liquids at temperatures at least a few degrees below the critical 
temperature. £ -> <» as T -* T c , the critical temperature. The two lines, GS, 

represent the Gibbs surface. 
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^4 ur ® 2 (i) Fourier optic* method of deeign for • eurface- 
light-scattering apectrometer. The symbol D implies a Fourier 
transform. Transmission (or reflection) functions relate 6 H to 
6***, s^ to s< 4 \ and d**’ to d* 41 . In particular, for tbs surface s C4> 
+ f + i»s w where v depends on the monolayer density 
and is ignorable. The phase function d depends on the ripple 
profile and is modeled easily. Ths paring transmits s zero-order 
beam of relatively high intensity which is passed along with higher 
order diffracted beams which can be filtered ipetially by Fl. (b) 
Schematic of an implementation for high-pressure work. Note 
that the grating along with filters F2 and Fl providt s heterodyne 
configuration. It la easy to block all beams at Fl except the 
aero-order beam so aa to obtain a bomodyne configuration. Note 
that the grating ia Imaged onto the surface. The polarising beam 
splitter (POLiUJ.) and A/4 plats Isolate the laser from the re- 
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